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摘 要 : 砾石 是 各 种 水 文 和 侵蚀 等 过 程 综合 作用 的 产物 ,是 草地 和 土壤 退化 ,生态 系统 恶化 的 一 个 标志 , 反 过 来 这 
些 砾石 也 影响 到 侵蚀 的 各 个 过 程 。 研 究 藏 北 高 原 地 表 砾 石 的 空间 分 异 对 区 域 生态 环 境 恢 复 具有 重要 意义 。 本 文 
以 地 表 砾 石 粒 径 大 小 和 空间 位 置 为 研究 对 象 ,通过 Moran’s 7 指数 .空间 变异 函数 .地理 探测 器 .回归 分 析 等 方法 对 
其 空间 异 质 性 进行 系统 性 分 析 。 结 果 表 明 :(1) 全 局 Moran’ s 1 指数 值 为 0.481, 旺 显著 的 正 相 关 , 局 部 Moran”s 148 
数 显 示 , 研 究 区 东部 砾石 聚集 模式 为 高 -高 聚集 ,中 部 为 低 - 低 聚集 ,其 余 区 域 多 呈 随 机 分 布 .(2) 砾石 空间 蜡 质 性 由 
结构 因素 主导 ,但 是 变异 函数 最 佳 拟 合 模型 与 特征 参数 值 均 存在 一 定 差异 性 , 即 存在 一 定 的 各 向 异性 特征 。(3) 地 


理 探 测 需 结果 显示 ,NDVI、 土 地 利用 类 型 为 影响 研究 区 砾石 粒 径 空间 异 质 性 的 主要 因素 ,人 口 密度 .植被 类 型 .年 均 
降水 为 次 要 因素 。(4) 回归 分 析 结 果 显 示 ,最 优 玉 度 回 归 为 最 佳 回 归 模 型 ,NDVI 对 砾石 粒 径 影响 最 大 ,其 后 依次 为 


土地 利用 类 型 .年 均 降 水 .植被 类 型 。 
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平均 海拔 4500 m 以 上 的 藏 北 高 原 ,是 我 国 的 重 
要 生态 安全 屏障 "。 但 是 ,近年 来 由 于 气候 变 暖 以 
及 人 类 活动 的 加 剧 , 藏 北 高 原 地 区 高 寒 草 地 退化 严 
重 , 植 被 覆盖 度 下 降 ,生态 环境 被 严重 破坏 ”, 这 直 
接 影响 着 藏 北 地 区 农 牧民 的 生产 生活 ,也 对 青藏 高 
原 乃 至 全 国生 态 造 成 严重 威胁 。 因 此 , 藏 北 高 原 
生态 环境 恢复 成 为 我 国生 态 工作 的 重要 内 容 ”。 而 
在 藏 北 高 原生 态 领 域 研 究 中 ,多 数学 者 以 草原 植被 
和 土壤 组 分 为 主要 研究 对 象 “" ,以 高 原 砾石 作为 
对 象 的 研究 相对 较 少 。 一 方面 ,砾石 是 各 种 水 文 和 
侵蚀 等 过 程 综 合作 用 的 产物 ,是 草地 和 土壤 退化 ， 
生态 系统 恶化 的 一 个 标志 ; 男 一 方面 砾石 反 过 来 又 
影响 水 文 和 侵蚀 过 程 ,如 入 渗 、 节 发 .径流 和 水 名 
等 中 。 因 此 ,本 文 以 地 表 砾 石 粒 径 大 小 和 空间 位 置 
为 研究 对 象 ,深入 分 析 砾 石 分 异 空间 格局 的 效应 关 
系 ,以 期 为 藏 北 高 原生 态 修复 提供 参考 。 

在 对 砾石 粒 径 空 间 分 异 的 研究 中 ,一 方面 在 研 
究 对 象 上 多 是 针对 沙丘 濒 海 ( 河 ) 的 沙 粒 以 及 戈壁 
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回归 分 析 


地 区 砾石 的 空间 分 异 , 如 Okin 等 在 加 利 福利 亚 东 
南部 的 马 尼 克 斯 (Manix) 研 究 了 羽 状 沙丘 沙 粒 空 间 
分 异 现象 ; 董 玉 祥 等 "对 河北 黄金 海岸 沙 粒 的 垂直 
分 布 模 式 进行 研究 ; 陈 西 庆 等 中 通过 对 长 江 入 河口 
底 沙 粒 径 的 变化 进行 研究 ,结果 表明 研究 区 河口 底 
沙 粒 径 存 在 继续 粗 化 的 趋势 。 另 一 方面 在 研究 方 
法 上 多 数 利用 经 典 统计 学 分 析 海 拔 .粒度 参数 与 粒 
径 的 相关 关系 ,或 利用 地 统计 学 中 的 空间 变异 函数 
模型 分 析 砾 石 粒 径 的 空间 异 质 性 ,如 Quade'5 通 过 
研究 位 于 葛 哈 维 戈壁 表面 的 砾石 分 布 随 海 拔 的 变 
化 ,构建 了 砾石 覆盖 度 随 海拔 变化 的 函数 ; 曹 晓 阳 
等 "9 利用 经 典 统 计 学 方法 分 析 了 噶 顺 戈壁 洪 积 遍 
区 地 表 砾 石 的 粒度 特征 ; 王 利 兵 等 加 采用 地 统计 学 
的 理论 和 方法 研究 了 浑 善 达 克 沙 地 沙 粒 粒 径 组 成 
的 空间 异 质 性 。 

上 述 研究 从 统计 学 及 空间 分 析 等 方面 对 砾石 
粒 径 空 间 分 异 展开 相应 研究 ,但 是 ,对 高 原 砾石 粒 
径 空 间 分 异 以 及 不 同 因子 春 加 贡献 率 的 研究 报道 
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较 少 。 地 理 探 测 融 模型 是 一 种 基于 空间 异 质 性 的 
研究 方法 ,可 以 定量 探测 表达 某 一 时 空 现象 的 主要 
驱动 因子 以 及 不 同 驱动 因子 之 间 交 互 作 用 ,与 其 他 
空间 异 质 性 探测 工具 相 比 地 理 探 测 吾 具有 更 高 的 
解释 效率 “。 因 此 ,本 文 在 借鉴 前 人 研 完 方 法 的 基 
础 上 ,加 入 地 理 探测 顺 模型 开展 藏 北 高 原 砾 石 粒 径 
空间 分 蜡 的 多 因子 定量 归 因 ,并 结合 回归 分 析 建 立 
砾石 粒 径 的 空间 预测 模型 ,以 期 为 藏 北 高 原生 态 修 
复 及 环境 保护 提供 科学 依据 。 


1 研究 区 与 研究 方法 


1.1 研究 区 概况 

由 于 藏 北 高 原 面积 较 大 , 且 部 分 地 区 砾石 分 布 不 
具 代 表 性 ,为 系统 分 析 藏 北 高 原 砾石 的 空间 异 质 性 ， 
故 研究 区 选择 人 砾石 分 布 特征 显著 的 那曲 地 区 南部 地 
X (29°55'21" ~ 33°53'22"N , 87°22'21" ~ 97°3'12"E) , 
包括 安 多 县 南部 ( 除 色 务 乡 外 其 余 安 多 县 境内 区 
域 )、 班 号 县 、 申 扎 县 、 双 湖 县 南部 ( 含 多 玛 乡 . 措 折 
罗 玛 镇 . 巴 岭 乡 、 协 德 乡 ) 以 及 色 尼 区 。 研 究 区 总 面 
积 约 18.31 km ,平均 海拔 4100 m DLE ,地势 四 周 高 、 
中 间 低 ,气候 寒冷 干燥 ,大 部 分 地 区 年 均 气 温 低 于 
0%。 年 均 降 水 量 从 东部 的 548.8 mm 左右 ,逐渐 减 


(a) 


。 地 级 市 m 研究 区 mam 那曲 地 区 mm 西藏 自治 区 


非 砾石 区 域 
(O) 


155km ™ 8~16 


至 西部 的 216.3 mm, 且 降 水 主要 集中 在 6 一 9 月 。 研 
究 区 土壤 以 寒 钙 土 和 草 秸 土 为 主 ,植被 类 型 较为 单 
一 ,90% 左 右 均 为 草 何 。 受 到 自然 地 理 条 件 的 影响 ， 
人 研究 区 经 济 发 展 水 平 较 差 ,经 济 发 展 以 放牧 业 与 旅 
游 业 为 主 , 人口 密 度 较 低 , 仅 在 东部 地 区 有 少量 人 
口 分 布 ,地 广 人 稀 。 
1.2 数据 来 源 

砾石 粒 径 空间 分 布 数据 来 自 野 外 采样 与 室内 
反 演 。 野 外 采样 点 位 于 藏 北 高 原 的 那曲 南部 地 区 
(图 1) ,海拔 每 升 高 100 m 设 一 条 样 带 , 每 个 样 带 设 
4~5 个 样 区 ,每 个 样 区 的 4 个 边 角 及 中 间 位 置 设置 1 
个 0.5 mx0.5 m 的 样 方 (图 1d), 由 于 区 域 情况 复杂 ， 
部 分 地 区 受 现实 条 件 的 限制 ,无 法 测 得 5 个 样 方 ,总 
计 测 得 30 个 样 区 (图 1b),145 个 样 方 。 每 个 样 方 内 
选取 占 主体 粒 级 10 个 以 内 的 砾石 ,并 用 游标 卡尺 对 
单个 砾石 的 Feret 直径 进行 测量 (图 1e) ,并 记录 海 
Th 坡度 . 坡 向 、 砾 石 数量 等 基本 信息 。 室 内 反 演 主 
要 将 145 个 样 方 分 为 130 个 反 演 样 方 和 15 个 验证 样 
方 , 通 过 采用 Image] 软件 计算 130 个 反 演 样 方 中 经 
形 貌 轮 廓 提取 后 的 野外 采样 砾石 的 粒 径 值 ,并 对 人 研 
究 区 Landsat 8 原始 影像 .光谱 指数 因子 .DEM 数据 
等 进行 了 预 处 理 并 去 除 因子 间 多 重 共 线 性 ,分析 粒 


EET, 


注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 CS(2019)3333 号 ,对 底 图 边界 无 修改 。 下 同 。 


图 1 野外 采样 与 反 演 成 果 


Fig. 1 Field sampling and inversion results 
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径 值 与 上 述 因子 间 的 相关 性 ,筛选 显著 变量 因子 ， 
结合 多 元 线性 逐步 回归 与 基于 PCA 的 多 元 回归 分 
析 构 建 砾石 特征 参数 反 演 模型 ,进而 得 到 砾石 粒 径 
反 演 数据 ”i。 该 数据 经 15 个 验证 样 方 验证 , 均 方 根 
误差 为 1.4, 相 对 误差 在 30% 以 内 的 样 区 数 为 11 个 ， 
约 占 总 样 区 数 的 73% ,平均 误差 为 22.75%, 反 演 精 
度 良 好 。 

根据 Wentworth 提 出 的 砾石 分 级 理论 ,将 本 次 
采集 的 砾石 粒 径 分 为 3 种 等 级 。 其 中 粒 径 在 2~4 
mm 之 间 的 砾石 分 为 极 细 砾 , 粒 径 在 4~8 mm 之 间 的 
分 为 细 砾 , 粒 径 在 8~16 mm 的 砾石 分 为 中 砾 。 本 研 
究 反 演 得 到 的 砾石 粒 径 分 布 如 图 1c 所 示 , 由 图 可 
知 ,研究 区 粒 径 主要 为 4~8 mm 的 细 砾 ,分 布 于 研究 
区 的 中 西部 , 占 总 面积 的 76.99% ,其 次 为 8~16 mm 
的 中 砾 ,主要 分 布 于 研究 区 东部 , 占 总 面积 的 
19.32% ,而 极 细 粒 仅 占 总 面积 的 3.69%。 

环境 变量 数据 有 自然 因素 、 人 文 因 素 两 大 类 ， 
共 8 个 代理 变量 ,数据 信息 如 表 1 所 示 。 
1.3 研究 方法 
1.3.1 数据 预 处 理 方 法 ”由 于 地 理 探测 器 的 正确 运 
行 需要 分 类 数据 , 故 对 土地 利用 类 型 .土壤 类 型 植 
被 类 型 采用 原始 分 类 ,其余 影响 因子 采用 自然 间断 
点 分 类 法 进行 离散 化 ,经 分 异 及 因子 探测 器 不 断 调 
试 ,最 终 确 定 NDVI\ 年 均 降 水 分 9 类 ,夜间 灯光 强度 
分 8 类 ,高 程 分 7 类 ,人 口 密度 分 4 类 。 

FHH ArcGIS 渔网 工具 将 研究 区 划分 为 10 kmx 
10 km 的 研究 单元 ,然后 采用 双 线 性 插值 法 提取 每 
个 研究 单元 内 砾石 粒 径 值 了 与 分 类 后 环境 变量 的 类 
型 量 X, 从 而 将 两 种 数据 的 属性 进行 匹配 ,最 后 将 提 
取 的 数据 带 入 地 理 探测 器 模型 进行 计算 涉 。 具 体 
流程 如 图 2 所 示 。 


1.3.2 空间 异 质 性 方法 

(1) 空间 局 域 异 质 性 Moran’ s 1 指数 可 以 分 析 
空间 局 域 异 质 性 , 即 分 析 某 区 域内 砾石 粒 径 与 周围 
砾石 粒 径 的 空间 自 相 关 程 度 。Moran 's7 指 数 介 于 
(-1,1) ,7 值 越 接近 于 1, 表示 研究 单元 正 相关 性 越 
强 ;1 值 越 接 近 于 -1, 表 示 研 究 单元 负 相 关 性 越 强 ;/ 
值 趋 近 于 0 时 ,表示 研究 单元 呈 随 机 分 布 *。 


n 


SS- DY- 
jo (1) 


PDA 


IP LT Moran’ s 1 指数;n 为 研究 砾石 单元 个 数 ;Y 
和 分别 为 第 i 个 砾石 单元 和 第 j 个 砾石 单元 的 粒 径 
值 ; 静 为 空间 权重 和 矩阵 ;8 为 观测 值 的 方差 ;7 为 粒 
EREE, 

(2) 空间 结构 异 质 性 为 进一步 分 析 空 间 自 相关 
的 范围 及 空间 结构 异 质 性 ,引入 空间 变异 函数 探究 
其 空间 结构 异 质 性 ,并 计算 各 向 同性 及 各 向 异性 下 
砾石 粒 径 的 最 佳 拟 合 模型 及 特征 参数 值 , 其 公式 为 : 

N(h) 
R(h)= aN -zz 站 (2) 

式 中 : 2(x,) 和 2(x,,,) 分 别 为 区 域 化 随机 变量 z 在 空间 
中 x 和 xi 处 的 值 ;N(h) 为 抽样 间隔 等 于 hh 时 的 点 对 
BX; RIA) AE PRA 

(3) 空间 分 层 异 质 性 砾石 的 空间 异 质 性 极 易 
受 环境 条 件 的 影响 ,为 此 本 文 引 入 地 理 探测 需 探 测 
影响 藏 北 高 原 砾石 粒 径 空 间 分 异 的 主导 因子 。 地 
理 探测 需 可 以 将 变量 进行 分 层 , 用 来 探究 空间 分 层 
异 质 性 ,并 揭示 变量 间 存 在 某 种 关联 的 空间 分 析 模 
型 1。 该 模型 被 广泛 应 用 于 土地 利用 、 区 域 经 济 、 
精准 扶贫 等 领域 ,是 一 种 新 型 空间 分 异性 分 析 工 


表 1 数据 来 源 及 介绍 


Tab.1 Data source and introduction 


数据 类 型 数据 名 称 数据 来 源 分 辨 率 

自然 因素 DEM(X,) 国家 地 球 系 统 科 学 数据 中 心 (http:www.geodata.cn) 30m 
年 NDVI(X) 中 国 科 学 院 资源 环境 科学 与 数据 中 心 (http:www.resdc.cn ) 1000 m 
土地 利用 类 型 (x) 司 家 地 球 系统 科学 数据 中 心 (http:www.geodatacn ) 30m 
土壤 类 型 (X,) 国家 地 球 系 统 科 学 数据 中 心 (http:www.geodata.cn) 1000 m 
植被 类 型 (万 ) 中 国 科 学 院 资源 环境 科学 与 数据 中 心 (http:www.resdc.cn) 1000 m 
年 均 降水 (X。) 中 国 科 学 院 资源 环境 科学 与 数据 中 心 (http:www.resdc.cn ) 500 m 

人 文 因素 AHRR) 中 国 科 学 院 资源 环境 科学 与 数据 中 心 (http:www.resdc.cn) 100m 
夜间 灯光 强度 (X;) 美国 国家 地 球 物理 数据 中 心 (https:www.ngdc.noaa.gov) 500m 
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2 数据 处 理 流程 
Fig. 2 Data processing flow 
具 ” 。 在 本 研究 中 用 以 探测 砾石 粒 径 的 空间 分 异 
以 及 环境 变量 因子 对 人 砾石 粒 径 的 解释 程度 ,用 g 值 
度量 ,表达 式 如 下 : 
XN, 
1- (3) 


St NAKA ACR o ARTE sh NEE 
的 分 区 ,h=1,2,3,…, 工 ;L 表 示 分 区 数目 。g 的 大 小 
反映 了 空间 分 异 的 程度 ,g 值 越 大 ,表示 空间 分 层 异 
质 性 越 强 , 反 之 则 空间 分 布 的 随机 性 越 强 。 


2 结果 与 分 析 


2.1 基于 Moran’s7 指 数 的 空间 局 域 异 质 性 

Moran’ s7 指 数 可 以 用 来 探究 研究 区 砾石 粒 径 
的 局 域 异 质 性 ,研究 区 全 局 Moran’ s 1 指数 为 0.481， 
2Z 值 高 于 2.5 且 通过 1% 的 置信 检验 ,表明 研究 区 砾 
石 粒 径 呈 显著 的 正 相 关 ,空间 聚集 性 较 强 ,存在 异 
质 性 。 局 部 Moran’ s 1 指数 结果 显示 (图 3) ,研究 区 
砾石 粒 径 聚集 模式 主要 为 高 -高 聚集 模式 、 低 - 低 聚 
集 模式 与 随机 模式 。 在 空间 分 布 上 ,研究 区 东部 砾 
石 多 呈现 高 -高 聚集 模式 ,研究 区 南部 砾石 则 以 随 
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/ æ 随机 分 布 
低 - 低 聚 集 m 高 - 低 聚 集 
O 研究 区 边界 = 水 体 m 高 -高 聚集 m 低 -高 聚集 
图 3 研究 区 砾石 粒 径 的 空间 聚集 特征 
Fig. 3 Spatial clustering characteristics of gravel size 


in the study area 


机 分 布 为 主 ,在 研究 区 西部 和 北部 地 区 ,随机 分 布 
和 低 - 低 聚集 呈现 相间 分 布 。 
2.2 基于 空间 变异 函数 的 空间 结构 异 质 性 

通过 空间 变异 函数 计算 出 研究 区 砾石 粒 径 的 
一 系列 变异 参数 ,并 选择 最 佳 空间 变异 函数 模型 进 
行 拟 合 , 从 而 得 到 研究 区 砾石 粒 径 的 空间 变异 曲线 
(图 4a)。 考 虑 到 砾石 粒 径 在 不 同方 向 上 的 变化 存 
在 一 定 差 异性 ,本文 从 各 向 异性 的 角度 出 发 ,将 研 
究 区 砾石 粒 径 划分 为 0" .22.5" 45° , 67.5° , 90° , 
112.5° 、135° 及 157.5° 共 8 个 方向 ,各 个 方向 上 的 点 
线 图 如 图 4b 所 示 。 

由 图 4 可 知 ,从 各 向 同性 的 角度 出 发 ,总 的 砾石 
粒 径 空间 变异 函数 类 型 为 线性 ,而 从 各 向 异性 的 角 
度 来 看 ,0° 、22.5° 45° .157.5" 方 向 上 变异 函数 类 型 


变异 函数 值 


图 4 砾石 粒 径 空间 变异 曲线 


Fig.4 Spatial variation curve of gravel size 
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为 指数 模型 90° 、112.5° .135° 方 向 上 变异 函数 类 型 
为 高 斯 模型 ,67.5° 方 向 上 为 线性 模型 。 

根据 变异 函数 曲线 ,得 到 研究 区 砾石 粒 径 空间 
结构 异 质 性 相关 特征 参数 如 表 2 所 示 。 

从 各 向 同性 角度 分 析 , 块 金 值 CG, 为 0.391, 具 有 
一 定 的 块 金 效 应 , 亦 说 明 随 机 因素 引起 的 空间 异 质 
性 不 可 忽略 。 偏 基 台 值 代表 结构 因素 引起 的 空间 
变异 ,为 0.708, 大 于 块 金 值 ,表明 研究 区 砾石 粒 径 的 
空间 异 质 性 主要 由 结构 因素 引起 。 基 台 值 为 1.099 ， 
代表 了 由 随机 因素 及 结构 因素 共同 引起 的 空间 变 
异 , 亦 表明 了 研究 区 砾石 粒 径 的 空间 异 质 性 不 可 忽 
略 。 块 基 比 为 0.356, 表明 研 究 区 砾石 粒 径 空间 异 质 
性 的 因素 中 ,随机 因素 占 比 0.356, 结 构 因 素 占 比 
0.644 ,结构 因素 占 主导 地 位 。 变 程 为 26.178 km, 表 
明 研究 区 砾石 粒 径 的 空间 自 相 关 范 围 为 26.178 kmo 

从 各 癌 异 性 的 角度 分 析 , 各 方向 特征 参数 值 存 
在 差异 性 ,但 在 任意 方向 上 , 偏 基 台 值 均 大 于 块 金 
值 ,说 明 研究 区 砾石 粒 径 空间 异 质 性 均 由 结构 因素 
主导 。 空 间 变 异 函 数 变 程 的 大 小 反映 砾石 粒 径 沿 
某 一 方向 变化 速度 的 快慢 , 变 程 值 越 大 ,表明 该 参 
数 在 该 方向 上 变化 慢 , 即 非 匀 质 性 越 弱 ,反之 越 
强 。 由 表 可 知 ,157.5" 方 向 变 程 值 最 大 , 即 该 方向 变 
异性 最 小 ,而 90° 方 向 变 程 值 最 小 , 即 该 方向 变异 性 
最 大 。 故 研究 区 砾石 粒 径 各 向 异性 特征 非常 明显 ， 
各 癌 异 性 主轴 为 157.5° 方 向 。 


表 2 砾石 粒 径 空 


2.3 基于 地 理 探 测 器 的 空间 分 层 异 质 性 

因子 探测 器 探测 各 因子 对 砾石 粒 径 的 解释 度 (0 
值 ) 大 小 。 探 测 结果 如 表 3 所 示 , 通 过 显著 性 检验 的 
因子 中 , 按 对 砾石 粒 径 解释 度 强 弱 排 序 为 :X(ND- 
V1)> 名 (土地 利用 类 型 )>( 人口 密 度 )>X;( 植 被 类 
型 )>Xs( 年 均 降 水 ), 对 应 的 g 值 分 别 为 0.41、0.27、 
0.16、0.15、0.13。 其 中 , NDVI 的 因子 解释 力 最 大 
(0.41) ,表明 藏 北 高 原 砾石 粒 径 空 间 分 布 受 NDVI 的 
控制 作用 最 为 强烈 , 即 NDVI 与 砾石 粒 径 在 空间 分 
布 上 具有 最 强 的 一 致 性 。 次 要 因素 为 土地 利用 类 
型 ,因子 解释 力 为 0.27 ,说明 土地 利用 类 型 也 是 影响 
砾石 粒 径 空 间 分 布 的 重要 因素 。 

交互 探测 主要 分 析 各 因子 对 砾石 粒 径 空间 分 
布 是 否 存在 交互 作用 。 结 果 如 表 4 所 示 ,由 表 可 知 ， 
除 成 站 和 表现 为 非 线 性 增强 外 ,其余 交互 作用 类 型 
均 表 现 为 双 因 子 增强 。 其 中 辫 (NDVDmR( 土 地 利 
用 类 型 ) 交 互 作 用 值 最 强 ,解释 力 达 到 0.56, 其 次 为 
X(NDVI)NX7 CA O i BE ) 5 X cE HHA FS NX. 
(年 均 降 水 ) ,解释 力 均 为 0.47, 最 小 的 为 总 (年 均 降 
JK OX, H ERE) ,解释 力 为 0.27。 各 个 因子 间 的 
交互 作用 共同 促成 了 研究 区 砾石 粒 径 的 空间 分 异 
格局 。 这 说 明 任意 两 个 环境 因子 交互 后 对 砾石 粒 
径 空间 分 布 的 因子 解释 力 均 会 显著 提升 ,也 就 是 说 
砾石 粒 径 空间 分 布 是 由 各 维度 环境 因子 的 共同 作 
用 形成 。 


间 变 异 函 数 特征 参数 


Tab. 2 Characteristic parameters of spatial variogram of gravel size 


方向 块 金 值 (Co) 偏 基 台 值 (C) 块 基 比 (CyCorC) 变 程 (a)/km 决定 系数 (RR ) 
各 向 同性 = 0.391 0.708 0.356 26.178 0.981 

0° 0.537 1.185 0.312 69.500 0.500 

22.5° 0.534 1.330 0.286 76.500 0.788 

45° 0.539 1.195 0.311 69.500 0.585 

67.5° 0.544 1.234 0.306 90.420 0.802 
各 向 异性 

90° 0.617 1.203 0.339 62.840 0.505 

112.5° 0.606 1.347 0.310 70.240 0.833 

135° 0.601 1.195 0.335 64.360 0.509 

157.5° 0.555 1.503 0.270 110.000 0.801 

RI 砾石 粒 径 因子 探测 结果 
Tab. 3 Gravel size factor detection results 
DEM(X,) = 4B NDVICX,) 土地 利用 类 型 (XK) 土壤 类 型 (和 %) 植被 类 型 (X;) FKA) AORE) 夜间 灯光 强度 (X) 


4 值 0.11 0.41 0.27 0.17 
P 值 0.35 0.00 0.00 0.99 


0.15 0.13 0.16 0.02 
0.00 0.00 0.00 1.00 
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RA 砾石 粒 径 影响 因子 间 交 互 作用 


Tab. 4 Interaction between impact factors of gravel size 


C A+B 结果 解释 
XX:=0.56 X,(0.41)+X;3(0.27) C >max( A,B 双 因 子 增强 
XmXs=0.46 X,(0.41)+X;(0.17) C >max( A,B 双 因子 增强 
XX.=0.46 X,(0.41)+X.(0.13) C >max( A,B 双 因 子 增强 
XmX=0.47 X,(0.41)+X,(0. 16) C >max( A,B 双 因子 增强 
XmXs=0.42 X,(0.27)+X;(0.17) C >max( A,B 双 因 子 增强 
XNX=0.47 X,(0.27)+X.(0.13) C>A+B 非 线性 增强 
XMNX=0.39 Xs(0.27)+X;(0.16) C >max(A,B 双 因 子 增强 
XNX=0.28 X;(0.17)+X<(0.13) C >max(A,B 双 因 子 增强 
XNX=0.26 X5(0.17)+X;(0.16) C >max(A,B 双 因 子 增强 
XNX-=0.27 X6(0.13)+X7(0.16) C >max( A,B 双 因 子 增强 


风险 探测 器 探测 各 因子 不 同 分 层 对 砾石 粒 径 
空间 分 布 是 否 存在 显著 性 差异 及 各 分 层 中 砾石 粒 
径 均 值 大 小 。 探 测 结果 如 图 5 所 示 ,NDVI 中 ,0.79~ 
0.92 分 区 中 粒 径 均值 最 大 ,0.1~0.18 分 区 粒 径 均 值 
最 小 ,总 体 粒 径 均 值 与 NDVI 呈 正 相 关 , 除 0.24~0.32 
与 0.61~0.79 与 0.79~0.92 分 层 外 ,其 余 区 间 均 存在 
显著 性 差异 。 在 土地 利用 类 型 中 , 玻 林地 中 的 粒 径 
均值 最 大 ,盐碱地 中 粒 径 均 值 最 小 , 绝 大 多 数 分 层 
间 均 存在 显著 性 差异 。 植 被 类 型 中 高 山 垫 状 植被 
中 的 粒 径 均值 最 大 , 除 高 山 垫 状 植被 分 层 外 ,其余 
均 不 存在 显著 性 差异 。 年 降水 量 中 ,年 降水 量 在 
470.4~548.8 mm 区 间 内 粒 径 均值 最 大 ,261.3~281.4 
mm 区 间 内 粒 径 均 值 最 小 ,总 体 两 者 呈正 相关 。 

人 口 密度 因子 中 ,人 口 密度 与 砾石 粒 径 两 者 呈 
现 一 定 的 正 相 关 。 但 是 ,NDVI 与 人 口 密度 均 呈 现 
“ 东 高 西 低 " 的 分 布 态势 ,两 者 具有 空间 分 布 上 的 相 
似 性 , 且 该 区 人 口 稀 少 , 人 类 活动 对 砾石 粒 径 的 影 
响 相 对 有 限 ,结合 因子 探测 器 探测 结果 ,NDVI 对 砾 
石 粒 径 的 解释 力 远 大 于 人 口 密度 ,NDVI 仍 为 影响 
砾石 粒 径 分 布 的 主导 因素 。NDVI 一 定 程度 上 代表 
了 植被 覆盖 度 ,NDVI 大 , 则 植被 覆盖 度 相 对 较 高 ， 
水 土 保持 能 力 强 ,流水 风力 等 外 力 侵蚀 作用 较 弱 ， 
则 砾石 粒 径 值 较 大 ;反之 ,砾石 粒 径 越 大 , 则 土壤 水 
分 的 入 渗 能 力 越 强 , 则 植被 覆盖 度 越 高 。 因 此 ,人 口 
分 布 与 砾石 粒 径 可 能 存在 某 种 假 相 关 , 即 两 者 空间 
分 布 的 相似 性 可 能 是 由 于 植被 等 外 部 环境 引起 的 。 
2.4 多 元 线性 回归 、 随 机 森林 回归 与 最 优 尺度 回归 
验证 

由 地 理 探测 器 的 因子 探测 器 及 交互 作用 探测 


注 : 字 母 A.B 分 别 代表 影响 因子 对 砾石 粒 径 的 解释 力 ,字母 C 代 表 AmB 对 砾石 粒 径 的 解释 力 。 


器 可 知 , 上 述 8 种 影响 因子 均 对 砾石 粒 径 有 不 同 程 
度 的 影响 ,但 是 DEM 土壤 类 型 .夜间 灯光 强度 P 值 
过 大 ,未 通过 0.05 的 显著 性 检验 , 故 将 NDVI. 土 地 
利用 、 植 被 类 型 .年 均 降水 、 人 口 密度 利用 多 元 线性 
回归 进行 拟 合 , 拟 合 结果 显示 各 因子 的 VIF 值 均 远 < 
5, 说 明 各 因子 间 不 存在 多 重 共 线性 ,模型 拟 合 调 整 
后 如 图 6a 所 示 , 尼 为 0.426。 将 上 述 8 个 影响 因子 加 
入 随机 森林 回归 模型 进行 拟 合 ,结果 如 图 6b 所 示 ， 
由 图 可 知 ,随机 和 森林 回归 模型 拟 合 尼 为 0.508 ,随机 
森林 的 尼 高 于 多 元 线性 回归 。 

考虑 到 影响 因素 中 有 土地 利用 类 型 植被 类 型 
等 名 义 变 量 , 利 用 多 元 线性 回归 与 随机 森林 回归 模 
型 进行 拟 合 ,无 法 准确 地 反映 分 类 变量 不 同 取 值 的 
距离 , 故 本 文 引 入 最 优 尺 度 回 归 进 行 拟 合 。 将 砾石 
粒 径 、NDVI.\ 年 均 降 水 、 人 口 密度 作为 数字 型 变量 ， 
土地 利用 、 植 被 类 型 作为 名 义 变量 输入 模型 , 结 
显示 最 优 尺 度 回 归 调 整 后 尼 为 0.491 , 拟 合 结果 好 于 
多 元 线性 回归 ,但 次 于 随机 森林 回归 , 且 人 口 密 度 
未 通过 0.05 的 显著 性 检验 ,原因 可 能 为 人 口 密度 与 
砾石 粒 径 存 在 某 种 假 相关 ,考虑 噜 除 人 口 密度 再 次 
进行 回归 ,回归 结果 如 表 5 所 示 。 

改进 后 模型 拟 合 精度 由 0.491 上 升 到 0.518, 通 
过 0.01 的 显著 性 检验 , 拟 合 结果 好 于 多 元 线性 回归 
与 随机 森林 回归 ,上 且 各 因子 均 通 过 0.05 的 显著 性 检 
验 ,但 该 模型 仅 用 于 分 析 研 究 区 的 砾石 粒 径 , 模 
型 为 : 

Y=0.578X, + 0.308X, +0.064X,+0.102X, (4) 

SUF YA WR BLAS BY tn HE SD ; Xo X Xs X a HA 
NDVI, 土地 利用 类 型 .植被 类 型 .年 均 降 水 的 标准 
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区 之 间 存 在 显著 差异 (P<0.05 )。 


图 5 影响 因子 分 区 粒 径 大 小 


Fig. 5 Influence factor partition particle size 


SY FE St FRU AE IE HH EK) IA O R BE 
有 较 显 著 的 水 平 ,通过 公式 (4) 系 数 绝对 值 大 小 可 
知 ,在 研究 区 影响 砾石 粒 径 大 小 的 因子 中 ,NDVI E 
地 利用 类 型 起 主要 作用 ,植被 类 型 .年 均 降 水 起 次 
要 作用 ,结果 与 地 理 探测 器 吻合 。 


3 讨论 


由 于 自然 环境 与 人 类 活动 的 影响 ,砾石 粒 径 存 


在 空间 异 质 性 。 本 研究 发 现 , 藏 北 高 原 砾石 粒 径 存在 
一 定 程度 的 空间 分 异 , 且 这 种 空间 分 异 是 由 结构 因素 
主导 。NDVI. 土 地 利用 类 型 是 影响 研究 区 砾石 粒 径 
空间 分 异 的 主要 因素 ,年 均 降水 .植被 类 型 人口 密 度 
是 影响 研究 区 砾石 粒 径 空间 分 异 的 次 要 因素 。 

地 理 探 测 器 探测 结果 显示 ,NDVI 与 土地 利用 
类 型 是 影响 藏 北 高 原 砾石 粒 径 空间 分 异 的 主要 
素 , 这 主要 是 由 当地 的 生态 环境 决定 的 (图 7)。 在 
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图 6 多 元 线性 回归 与 随机 森林 回归 拟 合 


Fig.6 Multiple linear regression and random forest regression fitting 


#5 最 优 尺 度 回 归 模 型 
Tab. 5 Optimal scale regression model 
改进 后 标 ”改进 后 
准 化 系数 ”显著 性 


标准 化 系数 ”显著 性 


NDVI 0.530 0.000 0.578 0.000 
土地 利用 类 型 0.313 0.000 0.308 0.000 
年 均 降 水 0.084 0.003 0.102 0.003 
植被 类 型 0.070 0.000 0.064 0.000 
人 口 密度 -0.044 0.212 一 S 


注 : 初 步 最 优 尺 度 回归 模型 调整 后 尼 为 0.491 ,改进 后 最 优 尺 度 回 
归 横 型 调整 后 尼 为 0.518。 


图 7 环境 变量 对 砾石 粒 径 的 影响 机 制 


Fig.7 The influence mechanism of environmental 


variables on gravel size 


研究 区 中 ,降水 .NDVI 呈现 “ 东 高 西 低 ” 的 态势 , 植 
被 类 型 土地 利用 类 型 以 高 寒 禾 草 、 苔 草草 原 ,高 寒 
斋 草 、 杂 草 类 草包 为 主 ,四 者 具有 空间 分 布 上 的 一 
致 性 。 降 水 影响 植被 覆盖 度 和 生长 状况 ,植被 类 型 
反映 了 植被 的 背景 信息 ,二 者 均 对 NDVI 与 土地 利 
用 类 型 产生 影响 。NDVI 反 映 了 植被 覆盖 度 与 植被 


生长 状态 ,其 与 土地 利用 类 型 又 与 水 流 、 风 力 等 外 
力 侵蚀 息息相关 (如 高 覆盖 度 草 地 与 低 覆 盖 度 草地 
相 比 ,水流 、 风 力 等 外 力 侵蚀 作用 较 弱 ) 。 而 砾石 一 
方面 是 各 种 水 文 和 侵蚀 等 过 程 综合 作用 的 产物 , 另 
一 方面 砾石 又 通过 入 活 改变 土壤 储 水 量 等 反 过 来 
影响 植被 和 土壤 肥力 ,进而 影响 NDVI 与 土地 利用 
类 型 。 

综 上 所 述 , 可 以 从 砾石 粒 径 角度 分 析 草 地 退化 
趋势 ,为 草原 牧草 退化 提供 动态 监测 。 人 类 活动 导 
致 藏 北 高 原 天 然 草 地 植被 急剧 下 降 ,中 高 覆盖 草地 
癌 低 覆盖 草地 转变 ,居民 和 工业 用 地 面积 增加 , 土 
地 利用 类 型 发 生变 化 ,导致 流域 水 土 保持 能 力 下 
降 ,土壤 破坏 更 为 严重 。 植 被 退化 与 土地 荒漠 化 的 
相互 作用 ,恶化 了 区 域 生态 环境 ,严重 浪费 了 区 域 
水 资源 。 由 于 供水 不 足 和 缺 水 ,灌溉 草地 的 牲畜 负 
重 过 重 , 草 地 退化 .荒漠 化 和 水 土 流 失 更 加 严重 , 形 
成 了 恶性 循环 。 上 述 变 化 将 会 使 水 土 保持 能 力 下 
降 , 外力 侵 蚀 日 益 严 重 ,砾石 数量 增多 , 粒 径 变 小 。 
但 受 研 究 资料 限制 ,缺乏 长 时 间 序 列 的 砾石 粒 径 分 
布 数据 ,因此 无 法 确定 砾石 粒 径 变化 状态 ,后 期 将 
在 时 间 尺 度 上 进一步 研究 砾石 粒 径 的 空间 分 异 ,以 
期 为 藏 北 高 原生 态 环境 修复 提供 科学 的 参考 依据 。 


4 结论 


(1) 研究 区 砾石 粒 径 呈 显著 的 正 相 关 ,空间 聚 
集 性 较 强 ,存在 异 质 性 ;在 研究 区 东部 砾石 多 呈现 
高 -高 聚集 ,研究 区 南部 砾石 则 为 随机 分 布 ,在 研究 
区 西部 和 北部 地 区 ,随机 分 布 和 低 - 低 聚集 呈现 相 
间 分 布 。 
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(2) 研究 区 砾石 粒 径 具有 明显 的 各 向 异性 特 化 研究 加 生态 学 报 , 2020, 40(3): 964-975. [Hao Aihua, Xue 
征 „H H 结 构 素 主导 > 从 变异 K 数 模 型 来 看 总 Peng et al. Study on gua and soil deg- 
的 变异 函数 模型 为 线性 其 余 方向 的 模型 为 指数 of ae grassland on the Qinghai-Tibet Plateau[J]. Acta 
ye pe et Ecologica Sinica, 2020, 40(3): 964-975. ] 
Er H 4| ë ZUN ¥ VLN 4 
i 斯 ` 线 HE , 从 变 程 来 看 ,157.5 za 站 FT 性 最 小 . [6] LiG Y, Mu J P, Liu Y Z, et al. Effect of microtopography on soil 
90" 方 向 变异 性 最 大 。 respiration in an alpine meadow of the Qinghai-Tibetan plateau{J]. 
(3) 因子 探测 器 显示 NDVI 土地 利用 类 型 为 主 Plant and Soil, 2017, 421: 10.1007/s11104-017-3448-x. 
影响 因素 ,植被 类 型 .年 均 降 水 与 人 口 密度 为 次 [7] IWE, EZ, 刘国华. 青藏 高 原 高 寒 草 地 地 上 植物 碳 积累 速率 对 
ois my HE AH <A: ae 生态 系统 多 功能 性 的 影响 机 制 加 . 植物 生态 学 报 , 2021, 45(5): 
要 影响 因素 ; 交互 探测 器 显示 各 因子 存在 明显 的 交 生态 系统 多 功能 性 的 影响 机 制 四 . 植物 生态 学 报 (5) 
ane JEFA RRK = 496-506. [Sun Jian, Wang Yi, Liu Guohua. The effect mechanism 
HH IR 效应 NDVIN 土地 利 用 类 型 的 交互 作用 Hx of carbon accumulation rate of above-ground plants on the multi- 
JEJ ,年 均 降 水 门人 Fae BE 的 交互 作用 最 小 5 风 念 探测 functionality of ecosystem in alpine grassland of Qinghai- Tibet 
右 表 明 不 同 分 级 区 间 各 因子 对 砾石 粒 径 空间 分 布 PlateaulJ]. Chinese Journal of Plant Ecology, 2021, 45(5): 496- 
具有 一 定 的 差异 性 ,其 中 土地 利用 类 型 具有 最 强 的 ever 
间 影 81 MF, POUT, 安 克 俭 , 等. 青藏 高 原 高 寒 章 原 与 章 旬 土壤 可 他 
区 间 影 响 差异 。 i 
Yr ane eerie PERU REF [J]. 水 土 保 持 通报 , 2021, 41(3): 87-93, 102. [Wei 
(4) 回归 分 析 结 果 显 示 ? ÈN E 回归 精度 最 Ning, Yu Wenzhu, An Kejian, et al. The key factors of soil erod- 
局， 人 口 密度 未 通过 显著 性 检验 ， 相关 系数 由 大 到 ibility of alpine grasslands and meadows on the Qinghai-Tibet Pla- 
小 为 :NDVI, 土地 利用 类 型 .年 均 降水 .植被 类 型 ; 且 teau[J]. Bulletin of Soil and Water Conservation, 2021, 41(3): 87- 
均 与 研 石 粒 径 旦 正 相关 , 除 年 均 降水 与 植被 类 型 与 。 31021 
地 理 探测 器 略 有 不 同 外 , 其 余 结 果 均 与 地 理 探测 器 [9] Wang Q F, Jin H J, Zhang T J, et al. Hydro-thermal processes and 
de BEL Sey ZS thermal offsets of peat soils in the active layer in an alpine perma- 
AAH y) a frost region, NE Qinghai-Tibet plateau[J]. Elsevier, 2017, 156(56): 
1-12. 
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Spatial heterogeneity of gravel size in Northern Tibetan Plateau 


XU Tao', YU Huan, KONG Bo’, QIU Xia, HU Mengke', LING Pengfei 
(1. School of Earth Sciences, Chengdu University of Technology, Chengdu 610059, Sichuan, China; 2. Chengdu 


Institute of Mountain Land and Disasters, Ministry of Water Resources, Chinese Academy of Sciences, Chengdu 
610041, Sichuan, China; 3. Sichuan Real Estate Registration Center, Chengdu 610014, Sichuan, China) 


Abstract: Gravel is the product of various hydrological and erosion processes and is a symbol of grassland and 
soil degradation and ecosystem deterioration. Consequently, gravel also affects various erosion processes. 
Studying the spatial differentiation of gravel in the Northern Tibetan Plateau is of considerable importance for the 
restoration of regional ecological environments. In this paper, the size and spatial location of surface gravel were 
studied, and the spatial heterogeneity was systematically analyzed by Moran’ s 7 index, spatial variogram, 
geographic detector, and regression analysis. The following results are presented. (1) The global Moran’ s Z index 
is 0.481, which shows a significant positive correlation. Meanwhile, the local Moran’ s Z index shows a high-high 
gravel accumulation pattern in the eastern part of the study area, low-low in the middle part, and mostly random 
distribution in the rest of the study area. (2) Gravel spatial heterogeneity is dominated by structural factors. 
However, some differences are found between the best fitting model of the variogram and the values of 
characteristic parameters; that is, certain anisotropy characteristics exist. (3) The geographic detector results 
showed that NDVI and land-use type were the main factors affecting the spatial heterogeneity of gravel size in 
the study area, while population density, vegetation type, and annual average precipitation were the secondary 
factors. (4) The results of regression analysis revealed that the optimal scale regression was the best regression 
model, and NDVI had the largest influence on gravel size, followed by land-use type, annual mean precipitation, 
and vegetation type. 

Keywords: Northern Tibetan Plateau; gravel size; spatial heterogeneity; geographic detector; regression 


analysis 


